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Abstract 

The type la supernovae (SNe la) observational data are used to estimate the parameters of a cosmolog- 
ical model with cold dark matter and the generalized Chaplygin gas model (GCGM). The GCGM depends 
essentially on five parameters: the Hubble constant, the parameter A related to the velocity of the sound, 
the equation of state parameter a, the curvature of the Universe and the fraction density of the general- 
ized Chaplygin gas (or the cold dark matter). The parameter a is allowed to take negative values and to 
be greater than 1. The Bayesian parameter estimation yields a = — 0.86lg '15 > Hq — 62.0tl'42 km/ Mpc.s, 
Q k0 = -1.26tJ;g, Q m0 = O.OOtg;^, Q c0 = 1.39±i;^, A = 1.00±°;°°, t = 15.3±f l and q = — 0.80±g;H, 
where to is the age of the Universe and qo is the value of the deceleration parameter today. Our results 
indicate that a Universe completely dominated by the generalized Chaplygin gas is favoured, what re- 
inforces the idea that the this gas may unify the description for dark matter and dark energy, at least 
as the SNe la data is concerned. A closed and accelerating Universe is also favoured. The traditional 
Chaplygin gas model (CGM), a — 1 is not ruled out, even if it does not give the best-fitting. Particular 
cases with four or three independent free parameters are also analysed. 



PACS number(s): 98.80.Bp, 98.80.Es, 04.60.Gw 



1 Introduction 



One of the most important problems today in cosmology is the nature of dark matter and dark energy that 
must dominate the matter content of the Universe £Q. The existence of dark matter is suggested by the 
anomalies in the dynamics of galaxies and clusters of galaxies Dark energy seems to be an inevitable 
consequence of the present acceleration of the Universe as indicated by the type la supernovae (SNe la) 
data, which asks for a fluid of negative pressure which does not agglomerate at small scales [310]. While in 
general the nature of dark matter can be connected with relics of fundamental theories, like axions [5], there 
are two main candidates to represent dark energy: cosmological constant 6 ; and quintessence which is a self 
interacting scalar field [Zj-IH]. Both models suffers of many drawbacks, mainly linked with fine tuning of 
parameters, as in the quintessence model |1(J) . or disagreement with the theoretically predicted values, like 
in the cosmological constant case 

The traditional Chaplygin gas model (CGM) and the generalized Chaplygin gas model (GCGM) have 
been widely considered as alternatives to the cosmological constant and to quintessence as the dark energy 
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that drives the present acceleration of the Universe ^S^dl- The appealing of the GCGM comes from, among 
other reasons, the fact that it can unify the description of dark matter and dark energy. The generalized 
Chaplygin gas (GCG) is characterized by the equation of state 



.4 



(1) 




where A and a are constants that in general are considered as positive numbers. The conservation equation 
leads to the relation 



Some of the main interesting features of eq. are the following. Initially, the GCGM behaves as a dust 
fluid, with p cx a~ 3 , while at late times the GCGM behaves as a cosmological constant term, p oc A 1 ^ 1+a \ 
Hence, the GCGM interpolates a matter dominated phase (where the formation of structure can happen) 
and a de Sitter phase. At the same time, the sound velocity of the GCGM is positive which assures that, 
in spite of exhibiting a negative pressure, it is stable against perturbations of the background configuration 



A lot of effort has been done in order to constrain the free parameters of the GCGM. Type la supernovae 
data ^7]-[2i>|, the spectra of anisotropy of cosmic microwave background radiation |27], the mass power 
spectrum j^Hl IH3 EH], gravitational lenses X-ray data [22] and age estimates of high-z objects J22J 
have extensively been used in this sense. This is essential for many obvious reasons. One of them is that 
if the parameter a is close to zero, the GCGM reduces essentially to the cosmological constant model. If 
observations indicate that this is the preferred value, the appealing for the cosmological constant as dark 
energy is reinforced. In a previous work |17|. we have analysed the case of the traditional Chaplygin gas 
model, where a = 1, using the SNe la data. One important aspect of this analysis is the use of the Bayesian 
statistics, for one side, and also the fact that all free parameters (Hubble constant, the curvature term, the 
value of the sound velocity and the proportion of ordinary matter and Chaplygin gas) of the model were 
taken into account in the treatment of the problem. This contrasts with most of the analyses previously 
found in the literature in the sense that they employ a more simplified statistical analysis or some input 
parameters were fixed a priori |18 | -|26 | . In that work |17j . it was concluded that the unification scenario 
(i.e., without pressureless matter), and a sound velocity near (but not equal) to the velocity of light were 
preferred. Moreover, a closed Universe was clearly favoured. 

The aim of the present work is to apply the same analysis, using the Bayesian statistics and taking into 
account all free parameters, to the GCGM. In order to keep contact with the previous work, we will use 
the same supernovae sample. This sample is much smaller than those available today (see, for example, ref. 
E3~|3ni)- But, the quality of the data is extremely good, which leads to a smaller value of xt (the quantity 
that measures the quality of the fitting): with this restricted sample the best-fitting gives xl around 0.74, 
while for the complete "gold" sample [21], the value mounts to unity. In the GCGM there are five free 
parameters instead of the four free parameters of the traditional Chaplygin gas model (a = 1): the Hubble 
constant Hq, the parameter A related to the sound velocity of the Chaplygin gas, the curvature density 
parameter O/-0, the ordinary matter parameter Q m o (or alternatively, the Chaplygin gas parameter f2 c o) and 
the equation of state parameter a. 

One important point in performing this generalized analysis is the range of validity of the parameter 
a. Usually, it has been assumed that Q < a < 1. The reason is twofold: a greater than one could lead, 
even if in the future, to a sound velocity greater than the velocity of the light; a negative should lead to an 
imaginary sound velocity, that is, instability. However, an analysis of the supernovae data may taken into 
account values of a greater than one or negative. In particular, in reference |25| . using a different statistical 
approach with respect to the method employed here, a value of a near 2 was favoured. In fact, values of a 
outside that limited range can be considered. The reason lies in the fact that the GCGM with the equation 
of state (pQ) may be a phenomenological manifestation of a theory expressed in terms, for example, of scalar 
field with a potential, or even of tachyon condensate model |37II38| . In this case, the notion of sound velocity 
can not be identified with the ordinary expression used in the perfect fluid approach, and no violation of 
causality or instability may occur if a is greater than one or negative. 




(2) 
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Hence, in our analysis the parameter a will vary from negative to highly positive values. The prediction 
for each parameter (a, Ho, flko, ^mo, ^co and A) is taken marginalizing, i.e., integrating over the other 
parameters. We stress that such marginalization procedure gives predictions that can be quite different from 
a two dimensional analysis, for example, where one searches for the simultaneous preferred values for any 
two parameters. This marginalization is in fact essential in order to have the correct prediction for a given 
free parameter. We will consider the following cases: all five parameters free; spatial curvature zero (four free 
parameters); pressureless matter given only by the barionic component (four free parameters); pressureless 
matter absent (four free parameters); spatial curvature zero and no pressureless matter component (three 
free parameters); spatial curvature zero and pressureless matter component given by the barionic component 
(three free parameteres) . For the more general case with five free parameters the predictions we obtain are the 
following: a = -0.86±fe, H = 62.0t\H km/Mpc.s, n k0 = -1.26±{H Ono = O.OQjftoo, «co = l-39±^, 
A = l.OO+Qgg. Moreover the predicted age of the Universe and the value of dc deceleration parameter today 
are to = 15. 3+32 and q Q = — 0.80^Q'g|, respectively. The results are consistent with, for example, the age of 
globular clusters |39j . Two crucial parameters are a and A, since if a = and/or A = 1 are favoured, the 
GCGM becomes equivalent to the cosmological constant model (ACDM). The value of the parameter A is in 
fact peaked at unity. But for reasons discussed later, this does not imply that this is the most favoured value; 
indeed it is just expected to be near unity. On the other hand, the parameter a is peaked in a negative value, 
but the spread is very large. In particular, the traditional Chaplygin gas with a = 1 can not be excluded 
or even disfavoured: in most of the cases this value of a is favoured with a probability around 50%. Hence, 
our analysis does not exclude the traditional Chaplygin gas model. 

This paper is organized as follows. In next section, the basic equations and quantities are exhibited. 
Section 3 explains the type la supernovae fitting. In section 4, we display the results of the Bayesian analysis 
of type la supernovae data with the GCGM. We present our conclusions in section 5. 



2 Description of the GCG model 

Regardless the possibility of unification of dark matter and dark energy using the GCGM, we will consider 
a model for the Universe today with two fluids: pressureless matter and the GCG. In doing so, our aim is 
to verify whether the unified model is the case favoured by observations. Moreover, even if in the GCGM 
there is no need of a dark matter component, there is baryonic matter in the Universe which is represented 
also by a pressureless matter. Hence, the dynamics of the Universe is driven by the Friedmann's equation 

a\ 2 k 8ttG ( \ 

-J +~ 2 = —[*» + *>) . (3) 
a 

p m + 3-p m = , (4) 
a 

fc+s s(*-£) = • (5) 

where p m and p c stand for the pressureless matter and generalized Chaplygin gas component. Dot means 
derivative with respect to the cosmic time t. As usual, k = 0, 1, — 1 indicates a flat, closed and open spatial 
section. In the conservation equation for the generalized Chaplygin gas component (JSJ) , the equation of state 
||TJ has been introduced. 

The equations expressing the conservation law for each fluid Q-© lead to 



PmO \ , . B 



1/(1+") 

r m - a 3 ' P c ={ A + a 3(l+a) | ' ( 6 ) 

Henceforth, we will parametrize the scale factor such that its value today is equal to unity, ao = 1. Hence, 
f ^ i/(i+") 

PmO and p c o =< A + B > are the pressureless matter and GCG densities today. Eliminating from the 
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last relation the parameter B, the GCG density at any time can be reexpressed as 



I- A 



i/(i+ Q ) 



Pc = Pco|A+- 5 ^ y | , (7) 

where A = A/p C Q. This parameter A is connected with the sound velocity for the GCG today by the relation 
v 2 s = aA. 

The luminosity distance is given by the expression 

D L = ^n , (8) 

a 

r*i being the co-moving coordinate of the source. Using the expression for the propagation of light 

ds2 = = dt2 -r^ > (9) 

and the Fricdmann's equation (|3J), we can reexpress the luminosity distance as 

D L = (l + z)S[f(z)} , (10) 

where 

S(x)=x (fc = 0), S(x)=sinx (k = 1), S(x) = sinhx (k = -1) , (11) 

and 

/(z) = i?oi {JW*' + 1)3 + n c0 [A + p + i)3(i+«)(i - A)]V(i+«) - n k0 (z> + rfyT* ' (12) 



with the definitions 

_8nGp m0 _8irGp c0 _ k 

"mo ^— 775- , " c o - — 5— 772 > > U^J 

d ±i 6 J± H 



where O m o + H e o + Qko = 1. The final equations have been also expressed in terms of the redshift z — — 1 + -. 
The age of the Universe today can be expressed as 



(14) 



to _ 

T J (1 + Z>){n m0 (z> + 1)3 + n c0 [A + (Z> + 1)3(1+«)(1 - A)]V(l+°0 - Q fe0 (z' + l)2}l/2 

where T = (100/i/o) x 10 10 , so to has units of years. 

The value of the decelerated parameter today, qo = — |§, is given by 

n m0 + n c0 (i - 3A) 

qo = g ■ ( 15 ) 

The value dj of the scale factor that shows the beginning of the recent accelerating phase of the Universe 
is given by the roots of the equation 

a = a{H + H 2 ) , (16) 
and is related to the redshift value Zi such that — = tt— or Zi — — 1 + 211 . 



3 Supernovae fitting 

We proceed by fitting the SNe la data using the GCG model described above. Essentially, we compute the 
quantity distance moduli, 

^ 5l0g Gd) +25 • (17) 
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Table 1: The best- fitting parameters, i.e., when \t is minimum, for each type of spatial section and matter 
content of the generalized Chaplygin gas model. Hq is given in km/Mpc.s, A in units of c, to in Gy and a, 
in units of qq. m and n are the constants of the relation log(l — A) = ma + n. 



and compare the same distance moduli as obtained from observations. The quality of the fitting is charac- 
terized by the \ 2 parameter of the least-squares statistic, as defined in Ref. pj], 



* -l^ a 2 , +(J 2 . ■ (18) 



ti i 15 * ne measurement error, a 2 nz i is the dispersion in the distance modulus due to the dispersion in 



In this expression, /Zq i is the measured value, /z i is the value calculated through the model described above, 
7t q i is the measurement error, <r^ 

galaxy redshift due to peculiar velocities. This quantity we will taken as 

d\ogD L 



Op: 



'-o z , (19) 



where, following Ref. |3lE0]. a z = 200 km/ s. 

In this article, we have used the same 26 type la supernovae in table 1 of ref. J7j, so we do not need to 
analyse and compare again the CGM with the ACDM (this work has been done in ref. J7j), but instead we 
can concentrate our efforts on the GCGM. Although nowadays there are approximately 200 SNe la available, 
the choice of these 26 are confirmed to be of excellent quality by analysing the low values of xt (the estimated 
errors for degree of freedom, i.e., \ 2 divided by 26, the number of SNe la), shown in tabled xt values range 
between 0.738 to 0.748, while other larger SNe la samples have been used in many cosmological models 
E2JES1^2, typically resulting xt near the unity. 

The best-fitting independent parameters and other dependent parameters (functions of the indepen- 
dent parameters) for the generalized Chaplygin gas model (GCGM) are given by table ^ for six different 
cases of spatial section and matter content. The first case (GCGM) has five free independent parameters 
(a, Hq, fimo, fico, A), the second case {GCGM : k = 0) has four free independent parameters (a, Ho,£l m o, A), 
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the third (GCGM : il m Q = 0) and fourth cases (GCGM : f2 m o = 0.04) have four free independent param- 
eters (a,H ,n o0 ,A), the fifth (GCGM : k = 0,_O m0 = 0) and sixth (GCGM : k = 0, f} m0 = 0.04) cases 
have three free independent parameters (a, Ho, A). 

It is important to emphasize that, for each case, all free independent parameters are considered simul- 
taneously to obtain the minimum of xt- So, for example assuming the GCGM with k = and £l m o = 0, 
if we ask for the best simultaneous values of (a, Hq,A) then the answer is given by last column of table 
^ However, in this example, asking for the best value of a by weighing (marginalizing or integrating) all 
possible values of (Hq, A) is another issue which is answered by the Bayesian estimations of the next section, 
not by best-fitting in n-dimensional parameter space whose results are listed in tabled 

It is remarkable that all independent and dependent parameters are physically acceptable : the cold dark 
matter density parameter range is 0.13 < £l m o < 0.17 (when it is not fixed a priori), the age of Universe 
today sits between 14.0 Gy and 14.4 Gy, the value of the scale factor that marks the beginning of the 
recent accelerating phase of the Universe goes from 0.76 oo to 0.80 ao, etc. 

But the values of a and A are somewhat a surprise : a is extremely large and A is as close to 1 as a is 
large. The minimization of xt was obtained in the following way : starting with an initial global minimum 
taken from the n-dimcnsional discrete parameter space (see next section), usually around (a, A) (3.0, 0.95), 
we increment a slowly at each step, the new step uses the last step parameter values as initial values to 
search the local minimum of xt (by using the function FindMinimum of the software Mathematica |42j) 
and so forth. So, for each value of a there is a value of A, the relation between A and a is simple indeed : 
A = 1 — io ma +™ or log(l — A) = ma + n, with m « —0.3 and n ps —0.25, or see the tableHfor each case of 
spatial section and matter content. 

The minima of xt as function of a change very slowly. For example, the first case (GCG) has the 
minimum of xt starting with 0.750 (when a = 3.5), decreasing to 0.744 (when a — 7) and reaching the 
best minimum 0.738 (when a = 42.04). Through all the evolution of a, the independent and dependent 
parameters have physically acceptable values, so low values of a are only barely worse than large values of 
a with respect to xt minimization. 

Slightly smaller values of xt favour the GCGM over the CGM (compare with table 3 of JIT]), independent 
of the spatial section type and matter content. The next section employs the Bayesian statistics to obtain a 
better comparison between these models and a more robust estimation of parameters. 

The Mathematica package used for the Bayesian analyses, developed by one of the authors (R. C. Jr.), 
is due to be publicly released. It comprises functions for calculation, marginalization, credible/confidence 
regions and intervals, maximization and visualization. 

4 Bayesian analyses of the cosmological parameters 

The same Bayesian statistics approach presented in section 3 of ref. ^2] was employed here to obtain the 
parameter estimations and answers for some hypothesis about the generalized Chaplygin gas model (GCGM) 
with up to five free parameters (a, Hq, Q m o, ^co, A), the age of the Universe to, the deceleration parameter qo 
and the value of the scale factor that shows the beginning of the recent accelerating phase of the Universe. 
The estimations have a central value where the one-dimensional PDF (Probability Distribution Function) is 
maximum and the positive and negative uncertainties are defined at a 2a (95.45%) credible (or confidence) 
level, see more details in section 3 of ref. ^7] . 

Each independent parameter estimation used an one-dimensional PDF obtained by marginalization, 
i.e., where (n — l)-dimensional integrals are computed for each value of the parameter in the case of a n- 
dimensional parameter space. For example, in the 5-parameters case (without any restriction on fi^o and 
f2 m o), the p(a \ /iq) is obtained through marginalization process described by 4-dimensional integrals of 
p(a, Hq, fl m o, fl c o, A | yito) over a 4-dimensional parameter space. 

An ideal calculation to compute the n-dimensional integrals would include an infinite number of samples 
of parameter space points with infinite volume, but in practical estimations we are limited to choose a finite 
region of the parameter space (such that outside it the probabilities are almost null) and a finite number of 
samples. With this criteria, we have chosen — 1 < a ^ 7 with resolution from 0.2 to 0.5. The lower resolution 
of 0.5 is employed in the regions where the PDF changes slowly with respect to a, which is interpolated 
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Fi gure 1: The graphics of the joint PDF as function of (a, Ho, A) for the generalized Chaplygin gas model when fixing k = 
and Qm,o = 0.04. It is a 3D density plot where the function (here the PDF) is rendered as semi-transparent colourful gas, 
or we can say that the plot is made by voxels (volume elements) analogous to pixels (picture elements). Where the PDF is 
minimum, the transparency is total and the colour is red, where the PDF is maximum then the voxel is opaque and the colour 
is violet; mid-range values are semi-transparent and coloured between red and violet (red, orange, yellow, green, blue, violet). 
The parameter ranges are —1 < a ^ 7, 55 ^ Ho $C 70 and ^ A ^ 1, with resolutions 0.2 X 1 X 0.05, respectively. The joint 
PDF of (a, Ho, A) for other cases (5 free parameters, 4 free parameters, etc.) show a similar 3D shape. 

to generate a higher resolution of 0.2 for two and one-dimensional PDFs. Analogously, 55 ^ Hq ^ 70 
with resolution of 1, ^ f2 m o ^ 1 with resolution of 0.05, ^ fi c o ^ 3.4 with resolution of 0.17 and 
^ A ^ 1 with resolution of 0.05. So there is a 5-dimensional discrete space of 23 x 16 x 21 x 21 x 21 
points (total of 3,408,048 points whose PDF values were calculated) as well as lower dimensional discrete 
spaces (fixing f2 m oi etc.), and the continuous PDFs in three, two and one-dimensional spaces are obtained 
by marginalizing interpolated functions built from the discrete points. We obviously discard the parameter 
space regions where the PDFs are not real numbers, as well as to is not a positive real number. 

Figure n illustrates the 3-dimensional PDF of (a, H ,A), using a 3D density plot composed by voxels 
(volume elements) with transparency and colour attributes, specific to the 3 free parameter case of the 
generalized Chaplygin gas model with k = and Sl m o = 0.04 (the other five cases with 3, 4 and 5 free 
parameters have similar 3-D shapes). It is clearly visible the non-Gaussian behaviour of the 3-dimensional 
PDF, with the orange and green 3D credible surfaces sharpening as a increases. Table H] tells us that the 
maximum is located at (a,H Q ,A) = (27.0,62.7,1 - 5.5 x 10~ 9 ) but it is located at a very narrow region 
(Hq is extremely peaked and 1 — 1 x 10~ 8 < A < 1 with non-negligible PDF) therefore it would be invisible 
even if the graphics had the range — 1 < a ^ 28 (because the needed resolution in the A axis would require 
approximately 10 8 voxels in this dimension and only one voxel would be opaque and violet among them). 
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Table 2: The estimated parameters for the generalized Chaplygin gas model (GCGM) and some specific 
cases of spatial section and matter content. We use the Bayesian analysis to obtain the peak of the one- 
dimensional marginal probability and the 2<r credible region for each parameter. Hq is given in km/Mpc.s, 
A in units of c, tg in Gy and Oj in units of do . 

Instead, in figure Q] we see a opaque violet region indicating a maximum near (a,H ,A) — (3.0,62,0.95) 
because there is non- negligible volume with high PDF values. The conclusion is simple: as a increases the 
high PDF regions narrow in width until becoming invisible, like the blue, green and yellow credible surfaces 
(corresponding to constant PDF levels) do in sequence as shown by figure ^ 

The above behaviour is important to understand the results of the marginalization processes. For ex- 
ample, when the marginalization is applied to figure ^ i- e ? resulting the bottom- right plot of figures [3 The 
integration of p(ce, Hq, A | fio) over the Hq parameter space yields p(a,A | /xo), the PDF maximum moves 
to lower values of a, (a, A) = (2.75, 0.95), because the high PDF regions for higher values of a have a small 
width with respect to Hq and weight less than larger regions due to the integration over Hq. Taking the 
marginalization process even further, i.e., marginalizing p(a, A \ fj,o) over the A space to produce p(a \ fj,o) 
illustrated by the bottom-right plot in figure the PDF maximum is once again at lower a, a = —0.30, 
because the high PDF regions for higher values of a have a small width now with respect to A and weight 
less than larger regions due to the integration over A. 

Table |21 is very comprehensive and resumes our results, listing all the Bayesian parameter estimations 
for six different cases of spatial section and matter content within the context of generalized Chaplygin gas 
model. The following sub-sections will analyse in detail this tabic and accompanying figures. More emphasis 
is spent into the analyses of the generalized Chaplygin gas parameter a, because it is the subject of constant 
debate. 

In figures and El the one-dimensional PDF shapes are matched against half-Gaussians to estimate the 
probabilities outside the left and right edges, in this way increasing the precision of the parameter estimations 
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without excessively enlarging the parameter space. 



4.1 a — the generalized Chaplygin gas parameter 

The parameter estimation of a is strongly related to the A parameter, i.e., the joint PDFs for the parameter 
space of (a, A) do not have a Gaussian shape, so the marginalization process changes the peak values and 
credible regions. Table [21 figures |21 and [3] are carefully analysed below with respect to a. In brief, the 
estimations for a are quite spread (i.e., they have large credible regions) and they differ a lot depending on 
if two or one-dimensional parameter space is used. 

4.1.1 Generalized Chaplygin gas model with 5 free parameters 

In figures the GCGM with 5 free parameters has a joint PDF peak of 0.948 at (a, A) = (-0.52,0.73) 
with 1 a (68, 27 %), 2 a (95, 45 %) and 3a (99, 73 %) credible region contours with PDF values of 0.717, 0.123 
and 0.002, respectively. The two-dimensional credible regions are not Gaussian-shaped, and their extremes 
give the estimations (a, A) = (-0.52>+^ 8 52 , 0.73±^) at la, (a, A) = (-0.52>+^ 8 52 , 0.73±g;^) at 2 a, and 
(a, A) = (-0.52> o^4 8 52 , 0. 73^0 57) at 3 a confidence levels. The symbol (>) in the upper credible value for a 
is due to the fact that the graphics is limited to — 1 < a ^ 7. 

The marginalization (integration) of p(a,A \ fj,o) over the A space yields p(a \ (j,q), shown in figure |3| 
where the credible intervals are a = — 0.85]toiJ at 1 a and a — — 0. 85^0^5 at 2a (the 3a PDF level cannot 
be seen with a ^ 7). The PDF peak of 0.383 at a = —0.85 is greater than the PDF value of 0.195 for a = 1 
(ordinary Chaplygin gas model), so the region with PDF greater than 0.195 has a CDF of 0.607, i.e., we can 
say that a = 1 is ruled out at a 60.7% (0.85 a) confidence level (or it is preferred at 39.3%). Table shows 
that positive values of a are more likely at a 64.1 % (0.92 a) confidence level. 

4.1.2 GCGM with 4 free parameters and k — 

The 4-parameter case with k = has a joint PDF peak of 1.165 at (a, A) = (2.09,0.96) with la, 2 a 
and 3a credible region contour levels of 0.695, 0.141 and 0.007, respectively. The two-dimensional credible 
regions give the estimations (a, A) = (2.09>+^ g 91 , 0.96±g;j£) at la, (a, A) = (2.09>+^ 9 91 , 0.96^°;^) at 2 a, 
and (a, A) = (2.09>^g 9 \ 0.96±g;^) at 3 a confidence levels. 

The estimation based on p(a \ /j,q) yields a — — I.OO^q'qq at la and a = — I.OO^qoo a t 2<t confidence 
levels. The PDF peak of 0.323 at a = -1.00 is greater than the PDF value of 0.205 for a = 1 such that 
a = 1 is ruled out at a 53.1 % (0.72 a) confidence level. Positive values of a are preferred at a 70.5 % (1.05 a) 
confidence level. 

4.1.3 GCGM with 4 free parameters and O, m o = 

The joint PDF has a peak value of 1.239 at (a, A) — (—0.51, 0.60) with la, 2 a and 3a credible region contour- 
levels of 0.378, 0.017 and 0.007, respectively. These levels indicate that the joint two-dimensional PDF is 
strongly peaked. The two-dimensional credible regions give the estimations (a, A) = (—0.51^0 45, o.6o±°:i2) 
at 1 a, (a, A) = (-0.51>j4 9 51 , 0.60^;^) at 2 a, and (a, A) = (-OiP^ 51 , 0.60± ;^) at 3 a confidence levels. 

The estimations based on one-dimensional PDF yield the same peak value (due to the highly peaked 
joint PDF), with a = — 0.51±g || at 1 a and a = — 0.51±g || at 2 a confidence levels. The PDF peak of 0.263 
at a = —0.51 is greater than the PDF value of 0.200 for a = 1, and we can say that a = 1 is preferred at a 
55.9% (0.77 a) confidence level. Positive values of a are preferred at a 76.2% (1.18 a) confidence level. 

4.1.4 GCGM with 4 free parameters and il nM = 0.04 

The joint PDF has a peak value of 1.257 at (a, A) — (—0.56, 0.60) with la, 2 a and 3a credible region contour- 
levels of 0.428, 0.020 and 0.005, respectively, which show a highly peaked behaviour. The two-dimensional 
credible regions give the estimations (a, A) = (-0.561^, 0.60±g;fg) at 1 a, (a, A) = (-0.56>^ 4 56 , 0.60tg;^) 
at 2 a, and (a, A) = {-OMtnlf 6 , 0.60^ ^) at 3 a confidence levels. 



9 



PDF of (a, A) PDF of (a, A) for k=0 




Figure 2: The graphics of the joint PDF as function of (a, A) for the generalized Chaplygin gas model. The joint PDF peak 
is shown by the large dot, the credible regions of 1 a (68, 27%) by the red dotted line, the 2 a (95, 45%) in blue dashed line and 
the 3 <r (99, 73%) in green dashed-dotted line. The cases for Q m o = are not shown here because they are similar to the ones 
with flmO = 0.04. 

The estimations based on one-dimensional PDF give a = — 0. 54^0 46 at lcr and a — — 0.541q'46 at 2 a 
confidence levels. The PDF peak of 0.280 at a = —0.54 is greater than the PDF value of 0.202 for a = 1, 
and we can say that a = 1 is preferred at a 50.9 % (0.69 a) confidence level. Positive values of a are preferred 
at a 73.6% (1.12 a) confidence level. 

4.1.5 GCGM with 3 free parameters and k — 0, fi m o = 

The joint PDF has a peak value of 1.336 at (a, A) — (3.11,0.95) with 1<t, 2a and 3<r credible region con- 
tour levels of 0.583, 0.115 and 0.007, respectively. The two-dimensional credible regions give the estimations 
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Figure 3: The PDF of a for the generalized Chaplygin gas model. The solid lines are the PDF, the la (68.27%) regions are 
delimited by red dotted lines and the 2a (95.45%) credible regions are given by blue dashed lines. The cases for Q m o = are 
not shown here because they are similar to the ones with Q,„o = 0.04. 



(a, A) = (3.11™, Q.95™) at 1 a, (a, A) = (3.11±|;I?, 0.95±°;gl) at 2 a, and (a, A) = (3.11>™, 0.95™) 
at 3 a confidence levels. 

The estimations based on one-dimensional PDF give a = 0.141^4 at 1 a and a = 0.14+^4 at 2 a 
confidence levels. The PDF peak of 0.222 at a = 0.14 is greater than the PDF value of 0.208 for a = 1, and 
we can say that a = 1 is preferred at a 63.2 % (0.90 a) confidence level. Positive values of a are preferred at 
a 78.8% (1.25 a) confidence level. 



4.1.6 GCGM with 3 free parameters and k = 0, Q. m0 = 0.04 

The joint PDF has a peak value of 1.367 at (a, A) = (2.75,0.95) with la, 2 a and 3cr credible region con- 
tour levels of 0.595, 0.113 and 0.007, respectively. The two-dimensional credible regions give the estimations 
(a, A) = (2.75™, 0.95™) at 1 a, (a, A) = (2.75™, 0.95™) at 2 a, and (a, A) = (2.75>™, 0.95™) 
at 3 a confidence levels. 

The estimations based on one-dimensional PDF give a = — O.SO^Q'yJ at 1 a and a — — 0.301q'yq at 2 a 
confidence levels. The PDF peak of 0.242 at a = -0.30 is greater than the PDF value of 0.214 for a = 1, 
and we can say that a = 1 is preferred at a 53.1 % (0.72 a) confidence level. Positive values of a are preferred 
at a 76.0% (1.18 a) confidence level. 

The non-Gaussian shape of the credible region contours in the (a, A) two-dimensional parameter space 
is the culprit of this difference for the estimations of a, so it is crucial to emphasize how the a parameter is 
estimated : from a joint two-dimensional PDF or an one-dimensional PDF (after marginalization of the joint 
two-dimension PDF). The most robust parameter estimation is based on the one-dimensional PDF, therefore 
giving small importance to regions in the joint two-dimensions which have a small joint PDF volume (CDF). 
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Figure 4: The one-dimensional PDF of A for the generalized Chaplygin gas model. The solid lines are the PDF, the la 
(68.27%) regions are delimited by red dotted lines, the 2a (95.45%) credible regions are given by blue dashed lines and the 
3a (99.73%) regions are delimited by green dashed-dotted lines. The cases for Q m o = are not shown here because they are 
similar to the ones with Q m o = 0.04. 



Comparing the estimations, the ones based on two-dimensional parameter space always give larger values 
for a, after marginalization the PDF is peaked at lower values of a. 

We also expect that a large sample of SNe la will narrow the credible regions of a, so the issue of a 
ordinary Chaplygin gas model (CGM) being preferred or ruled out will be probably better addressed. 

4.2 The parameter A 

The generalized Chaplygin gas model (GCGM) estimations for A are given by table El which if compared 
with the results of the CGM (table 5 of ref. ^3) clearly indicate that the PDF peak is closer to A — 1 and 
the dispersion is increased when the GCGM is considered. This behaviour is due to the shape of the joint 
PDF in the parameter space (a, A), figures [21 explicitly show that as a increases there is a narrower region 
of high PDF values closer to higher values of A. 

But the estimations for A shown in table and figures have to be interpreted. They do not indicate 
that the ACDM limit (A = 1) is favoured because the resolution for the A sampling is 0.05, i.e, a PDF peak 
at A = 1.00 means the peak happens for 0.95 < A ^ 1 indeed. Therefore, in figures 01 the graphics for the 
PDF of A which are peaked at A = 1 have the shape of the graphics for the case £l m o = 0.04, with the PDF 
peak closer to A = 1. See ref. for a detailed comparison between the CGM and ACDM cases. 

The analyses of a presented above have considered the joint PDFs for the two-dimensional parameter 
space of (a, A), and the extremes of the credible regions have given estimations for (a, A) at la, 2a, and 
3 a confidence levels. These estimations for A are almost always different from the estimations based on the 
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one-dimensional PDF of A, obtained by marginalizing p(a, A \ /iq) over the a space. We stress again that 
the most robust estimations are based on the one-dimensional PDFs. The general behaviour of table El is 
that the PDF peak is at 0.95 < A ^ 1 and the 2 a credible interval begins at A ~ 0.5 and finishes at A ~ 1. 

4.3 The Hubble parameter H 

The parameter estimation of Hq shown in tablets almost the same for all cases, around 62 km/Mpc.s with 
a narrow credible region of ±3 km/Af>c.s. This estimation, like the one for the CGM in table 5 of ref. |17j . 
is compatible with other SNe la analyses. 

4.4 The curvature density parameter D,ko 

Tableland figure[5]show that the curvature density parameter O^o estimation for the GCGM is almost like 
for the CGM of ref. J7|, i.e., it also strongly depends on the cold dark matter content of the GCGM. 

In the 5-parameter case a closed Universe (fifco < 0) is preferred at 97.19% (2.20 a) and a spatially flat 
Universe is ruled out at 93.44% (1.84 a) confidence level, i.e., the PDF of flko = is smaller than the PDF 
of -2.57 < flko < and this region has a CDF of 93.44% (1.84tr). We can also see the figure that the 
integral of the joint PDF over the upper region gives the same high likelihood of a closed Universe, although 
the joint PDF peak is located in the open Universe region. 

But assuming the hypothesis that f2 m0 = 0.04, a 4-parameter case, the behaviour is opposite, a closed 
Universe (fifco < 0) is preferred at 57.75% (0.80 a) and a spatially flat Universe is preferred at 90.32% 
(1.66 a) confidence level, i.e., the PDF of £lko = is smaller than the PDF of < fi^o < 0.17 and this region 
has a small CDF of 9.68%. Therefore, if we impose a low value for r2 m o then the spatially flat Universe is 
favoured. 

4.5 The cold dark matter density parameter Q m0 

The generalized Chaplygin gas model favours small values of Q m o like the ordinary CGM, see tableland 
table 5 of ref. ^3- The graphics for O m o in figure for the 5-parameter case and 4-parameter case with 
k = 0, shows a estimation for f2 m o peaked at at a high level of confidence, with f2 m o < 0.42 (for the 
5-parameter case) and £l m o < 0.20 (for 4-parameter case with k = 0) at 1 a (68.27%) confidence level. As 
f2 m o was sampled with bin size 0.05, then the peaks can be somewhere between 0.00 and 0.05. 

The joint PDF for the parameter space of (r2 m o,r2 c o)j i-e, p(f2 m o,^ c o I Mo) ; is shown in figure 
The PDF peak is now located at (n m o,tt c o) = (0.0,0.81), with the same the central value of O m o af- 
ter Bayesian marginalization. But, in this two-parameter space with simultaneously values of (Jl m o,fiA), 
the credible regions would be given by (O m0 ,f2 c o) = (O.O^qq, 0.81+oVg*) at la (68,27%) and (fi m0 ,^co) = 
(0.0l 0.81±g;^) at 2 a (95, 45%) confidence levels. Of course, the integral (marginalization) of p(Q m o, fi c o | 
/jLo) over the f2 c o space yields p(fl m o | /io), shown in figure The conclusion is that the marginalization 
decreases the confidence regions, even if it does not change the peak value of probability. 

4.6 The generalized Chaplygin gas density parameter Q c0 

The generalized Chaplygin gas density parameter f2 c o estimation (table [5J is like the estimation in the case 
of CGM shown by table 5 of ref. ^7], and is quite spread, i.e., it has a large credible region, for example 
fl C Q = 1-391^25 m t ne case of 5-parameter space. Figure |B] show the PDF of Q c q for the 5-parameter 
space (showing n cQ = 1.39±gJ| at 1 a confidence level) and the 4-parameter space with fc = (when 
fl C Q = 1 — Sl m o). 

As discussed above for the GCGM, the joint PDF for the parameter space of (fi m o,r2 c o) in figure 
shows that the credible regions would be given by {n m0 ,n c0 ) = (O.Ot8:$,0.8ltJ:g) at \ a (68,27%) and 
(n m0 ,n c0 ) - (0.0tJ;^,0.8llj^) at 2(7 (95,45%) confidence level. So, the marginalization of p(Q m Q , Q c q \ 
/io) over the Q m o space gives p(fl c o \ ^o), yielding a different peak value as well as different (narrower) 
confidence regions. 
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Figure 5: The graphics of the joint PDF as function of (Q m o, £2 c o) for the generalized Chaplygin gas model, where p(Q m o, £l c o) 
is a integral of p(ct, Ho, ilmO, f2 c o, A) over the (a, Ho, A) parameter space. It shows a different shape for the confidence regions 
with respect to the ACDM but similar to the CGM of ref. 1171 . The joint PDF peak has the value 1.39 for (f2 m o,^co) = 
(0.00, 0.81) (shown by the large dot), the credible regions of 1 <r (68, 27%, shown in red dotted line), 2 a (95, 45%, in blue dashed 
line) and 3 a (99, 73%, in green dashed-dotted line) have PDF levels of 0.50, 0.13 and 0.02, respectively. As Q^o + ^m0 + ^ic0 = b 
the probability for a spatially fiat Universe is on the line f2 m o + t2 c o = b above it we have the region for a closed Universe 
(k > 0, Qfro < 0)i an d below, the region for an open Universe (k < 0, Q^o > 0). 

4.7 The present age t of the Universe 

The estimates for the dynamical age of the Universe today, to, given in tabled show that the generalized 
Chaplygin gas model increases the 2 a credible regions (e.g., to — 15.3^2 Gy), when compared to the 
estimations of the CGM in table 5 of ref. ^7] (e.g., to = 14.2if 5 Gy). The peak values of to range between 
14.2 Gy to 15.3 Gy, while the CGM features to peaked between 13.1 Gy to 14.8 Gy. We expect that a large 
sample of SNe la will narrow the credible regions and allow, by means of independent estimations of to from 
different observations (globular agglomerates, etc), to verify which cosmological model is preferred : GCGM, 
CGM and ACDM. 
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Fi gure 6i The one-dimensional PDF of f2fco, QmO and n c o for the generalized Chaplygin gas model. The solid lines are the 
PDF, the la (68.27%) regions are delimited by red dotted lines, the 2a (95.45%) credible regions are given by blue dashed 
lines and the 3a (99.73%) regions are delimited by green dashed-dotted lines. As Q c o = 1 — Q^q — OfeOi f° r ^mO = we have 
Q c q = 1 — Slfco, for O m o = 0.04 then Q c o = 0.96 — QkO an d for Q^o = we also have O c o = 1 — fimO- 

Tabled shows the to for the best-fitting parameters without any margmalization (i.e., the n-parameter 
model has its n parameters taken simultaneously yielding a minimum x 2 ) : the to ranges from lA.OGy to 
14.4 Gy, well inside the estimated credible regions for to- 



4.8 The deceleration parameter q of the Universe 

All the estimations based on the generalized Chaplygin gas model (table |2J suggest an accelerating Universe 
today, i.e., qo < 0, with a high confidence level ranging from 1.66 a (90.24%) to 3.72 a (99.98%), depending 
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on the curvature and cold dark matter densities of the Universe. The peak values of go range between —0.89 
to —0.59. The data of CGM in table 5 of ref. {TJ. show a similar behaviour. 

Although the go values of Tabled from —1.00 to —0.80, are inside the estimated credible regions for go, 
the estimated peak values of go differs from the go for the best-fitting parameters. 

4.9 The inflexion point a, for the scale factor 

We have also estimated the value of the scale factor that shows the beginning of the recent accelerating 
phase of the Universe. Table shows that the most probable value for a% ranges between 0.57 ao to 0.77 do, 
where ao is the value of scale factor today, with usually an small upper credible region and a larger lower 
credible region. 

The probability of a.; < ao is theoretically the same having an accelerating Universe today (go < 0), 
and these independent probability calculations agree almost exactly as table [5] shows, indicating that the 
Bayesian probability analyses of this work are accurate and reliable for the given SNe la data. 

We can obviously express the scale factor value in terms of redshift value a^/ao = 1/(1 + z^) or 
Zi = — 1 + (ao/a,i). For example, a; = 0. 67^037 ao corresponds to Zi = 0.49^q4^. We expect that a large 
sample of SNe la with large redshift values will help improve this type of estimation of Zi or equivalently a^. 

The a, values are almost equal in table ^for the best-fitting parameters, between 0.756 ao and 0.792 ao, 
and this range fits well inside the the estimated credible regions for a^ of table [21 

4.10 The eternally expanding Universe 

The probability of having a eternally expanding Universe, given here by demanding d > during the future 
evolution of the Universe, is estimated to be exactly one or almost one (at almost a 6a level), see table |5J 
The estimations for the CGM and ACDM cases (tables 4 and 5, ref. ^7]) were miscalculated, they are also 
exactly one or almost one. 

5 Conclusion 

In this paper we investigated some constraints on the generalized Chaplygin gas model (GCGM) using type 
la supernovae data and Bayesian statistics. The model is built with the Chaplygin gas and a pressureless 
fluid as the matter content of the Universe, whose density parameters today are f2 c o and O m o respectively, 
such that the curvature term is £lko = 1 — ^mo — ^co- The parameter A related to the sound velocity of 
the Chaplygin gas, the Hubble parameter Ho and the a parameter are also independent free variables. We 
allowed the parameter a to vary between negative and large positive values. In principle, this could lead to 
instability problems at perturbative level (a negative), or superluminal sound velocity (a greater than one), 
but both potential problems can be resolved using a fundamental description for the GCG. 

When the five independent parameters are free, the Bayesian estimation yields a = — 0.861q°5, Ho = 
62.0tiilkm/Mpc.s, n k0 = -1.26±l;g, O m0 = O.OOt&jg, n c0 = 1.3911$, A = 1.001°;°°, t = I5.3±f 2 2 , 
go = — 0.80to;o2 an d a i — ^-^-0^7 a o where t is the age of the Universe, go is the value of the deceleration 
parameter today and a^ is the value of the scale factor that marks the start of the recent accelerating 
phase of the Universe. We have also estimated the cosmological parameters fixing the pressureless matter 
(through the constraints of nucleosynthesis) or the curvature (imposing the flat condition) , or even both. The 
estimations of the remaining free parameters do not change dramatically in each case, with some exceptions. 
For example, the prediction for the curvature parameter depends on whenever the pressureless matter is 
fixed or not : a spatially flat Universe is favoured only if f2 m o is fixed at low values. 

In particular, we verify that the a parameter is strongly related to the A parameter. In this manner, 
the best value of a have the peak values and the credible regions changed in function of the marginalization 
process used. The a estimations are very spread and differ a lot if one or two-dimensional parameter space 
is used. Many authors ^B] _ |2E1) prefer to estimate a and A simultaneously using two-dimensional credible 
regions and they do not marginalize the PDF to obtain the more statistically robust one-dimensional PDFs for 
a and A, so any comparison should be taken carefully as the credible region estimations on two-dimensional 
parameter space are usually different from the one-dimensional estimations. 
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Here we take the marginalization approach to obtain one-dimensional PDFs, and one important point 
arises in this analysis : despite the dependency on the A parameter, the value a = 1, that represents the 
ordinary CGM, is not ruled out since this particular case is preferred at 39.3% to 63.2%, depending on the 
spatial section and matter content. 

We calculated the value of the scale factor when the dark energy starts to prevail, m = 0.67 ao, for the 
case of five free parameters. In general, the value for a* ranges between 0.57 ao and 0.77 ao, where ao is the 
value of the scale factor today. This results are coherent with the suposition that the recent acceleration 
of the Universe could be started with the redshift parameter z between and 1 obtained by the type la 
supernovae data. 

In the case of the ordinary matter parameter fi m o our analyses show that the preferred values are small, 
around zero. This situation is the same of the case of the ordinary Chaplygin gas. The result reinforces the 
idea of the Quartessence model that unifies dark matter and dark energy into a single fluid, the 

Chaplygin gas. The predicted age of the Universe is well above the age of globular clusters. At same time, 
the deceleration parameter is consistent with the estimations from the WMAP data. 

Further developments of the present work include the use of a large sample of type la supernovae. 
However, the sample of 26 SNe la used here has an excellent quality with respect to the low values of xt- 
Moreover, the dispersion of the parameter estimations obtained here has the same order of estimations using 
larger samples. In some cases our estimation are even less spread, as our preliminary investigation on larger 
SNe la indicates. 

We also plan to apply the Bayesian analysis to the GCGM using large scale structures data (2dFRGS) 
and to the anisotropy of the cosmic microwave background. In this way, we may consistently cross the 
estimations of different observational data and constrain more sharply the predictions to each free parameter 
of the model. 
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